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Ph| Abstract 
> ■ 

^ I We investigate the total effect of correlations on photoionization of atomic states with 

^ ' nonzero orbital momentum, in the nonrelativistic high energy asymptotic limit, consid- 

^ O^ . ering the exclusive case of the dominant final state of an initial neutral atom. We find 

I that the substantial cancellation of the dominant intra-shell correlations, which had been 

■ reported earlier, can be understood utilizing the closure properties satisfied by the eigen- 

■ functions of the nonrelativistic Hamiltonian. Considering the sum of correlations with all 
OS . states, occupied or not, we show that complete sum is equal to the contribution of the 

I high energy part of the continuum. Consequently there is a total cancellation between 

• ■ the contributions of the bound states and the low energy part of the continuum states. 

. This means that the real correlations in the physical atom, due to the sum rule over the 

Q\ \ occupied states, can be also obtained as the negative of the total contribution of low en- 

■ ergy bound and continuum unoccupied states. We calculate this in the framework of the 
J> . quantum defect model. As we would expect, the results are close to those obtained earlier 

^ I in particular cases by direct summation over the occupied states. However this approach 

' also allows us to see that the dominant intrashell correlation is going to be cancelled. 

We can also obtain some limits on the correlation effects by considering calculations with 
the screened Coulomb functions. The role of correlations in the exclusive photoionization 
processes is discussed, also the modification of correlations in the case of atomic ions. 



1 Introduction 

In this paper we study the total effect of the final state correlations on the amplitudes and 
cross sections of photoionization of atomic states with nonzero values of orbital momenta in 
the nonrelativistic high energy asymptotic limit. It is known that in this situation correlation 
effects still persist even in the high energy limit. We show that the sum of such correlations 
shows a tendency to cancel. The experimental data for photoionization of p states of external 
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shells of neon and argon by photons of the energies of about 1 keV [1], [2] can not be interpreted 
in the framework of the Independent Particle Approximation (IPA). In IPA the p electron is 
ionized by direct interaction with the photon. The authors of [1] suggested a mechanism of IPA 
breaking by the final state electron interactions. The photon ionizes rather the s electron of the 
same subshell. In a next step the outgoing electron moves the p electron to fill the hole in the s 
state created by the photon. The final state correlations have been studied later in [3]- [5], with 
the initial state IPA breaking effects being included in [4], where they were shown to be small 
except some special cases. It was understood that at still larger energies, much greater than all 
the binding energies of the atom, there is a large cancellation between correlations with various 
shells. Such cancellations were first found in the angular distributions of photoionization [3], 
and later at the amplitude level [5]. 

Since calculations involving cancellations require more precise knowledge of wave functions 
for the description of the bound electrons, we try to demonstrate the cancellations in another 
way. We consider asymptotic energies of the outgoing electron, 

S > /, 

where / is the ionization energy of the single-particle ground state, and we seek to obtain the 
asymptotic amplitude. Thus we assume E to be much larger then all single-particle bound 
state energies. Our analysis is completely nonrelativistic, i.e, we assume E <^ m with m for 
the electron mass (we employ the system of units with h — c = 1). We consider only relatively 
light atoms, with not very large values of the nuclear charge Z, describing the bound electrons 
by nonrelativistic functions, with corrections of the order {aZy being neglected. We focus on 
the case of p electron photoionization. 

Assuming that all initial electrons arc moving in the same self-consistent field, we show 
that these cancellations can be understood utilizing the closure properties satisfied by the 
eigenf unctions of the nonrelativistic Hamiltonian. Considering the sum of correlations with all 
states, occupied or not, we demonstrate that complete sum is equal to the contribution of the 
high energy part of the continuum. Thus there is total cancellation between the contributions 
of the bound states and the low energy part of the continuum states, for which we will give a 
precise definition. Hence the sum of the correlations can be expressed as the negative of the 
sum over low energy unoccupied states. 

We perform some explicit calculations for real atoms by calculating the contribution of 
low energy unoccupied states, using the quantum defect model combined with the Fermi- 
Segre theorem, and making a rough estimate of the contribution of the low energy part of the 
continuum. In our approach all such terms are positive. Our results are close to those obtained 
by direct summation over the occupied shells. These direct terms occur with both signs, in the 
cases where such summations were carried out [3, 5]. 

We have carried out direct calculations here for the Coulomb case, and we find certain limits 
on the correlation effects in this case. Using perturbative treatment of the screening we show 
that the magnitude of cancellations in the real atom is greater. 

In our analysis we have used the perturbative approach to the final state interactions of 
the electrons developed in [6] . This approach was employed earlier for investigation of the IPA 
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breaking in photoionization [4]. Inclusion of these effects removed or strongly diminished the 
discrepancy between experimental data and the IPA calculations. We shall use the expressions 
obtained in [4] throughout the paper. 

In particular calculations we use the simplifying assumption that overlap matrix elements 
between the orbitals of different subshells in the initial state ion and the final state ion is 
small and can be neglected. In this approximation the inclusive cross section coincides with 
the exclusive one in which the state of the spectator electrons does not change. The latter 
cross sections correspond to the experiments [1, 7]. We show that the correlations considered 
in the paper can also manifest themselves stronger in inclusive processes of photoionization 
accompanied by excitation of external electrons. In the case of atomic ions there will be less 
cancellation among the correlations. 

The paper is organized as follows. In Sec 2. we recall the main equations for the perturbative 
treatment of the IPA breaking effects in photoionization. In Sec. 3 we write the sum rules 
provided by closure and show that they have the consequence that there is total cancellation 
between the sum over all the bound states and the low energy continuum states. In Sec. 4 we 
obtain the correlations of the occupied states, both directly (from previous work [5]), and as the 
negative of the sum over unoccupied low energy states. In Sec. 5 we make explicit calculations, 
using some simplified models. We show the results of these approaches in Table 2, including 
also the cases of large Z for illustrative purpose. For the direct calculations we investigate the 
Coulomb case in Sec. 6, with the results given in Table 3, and give some discussion of screening 
in an effective charge approach. In Sec. 7. we consider the role of correlation in exclusive 
processes, and for atomic ions. We summarize in Sec. 8. Some details of computations are 
presented in Appendices. 



2 Perturbative treatment of IPA breaking effects 

We recall the general points of our perturbative approach [6] , restricting ourself here to asymp- 
totic analysis. We shall use a simplifying assumption that the bound electrons are described 
by single-particle wave functions. This is not a necessary assumption, and the approach has 
been employed for the case of correlated functions as well [8], [9]. Consider the asymptotic am- 
plitudes Fi for ionization of initial state i with quantum numbers i — n,i,iz- The final state 
interactions between the outgoing electron and the residual ion in their lowest order in the the 
amplitude beyond the independent particle approximation (IPA) can be expressed in terms of a 
linear combination of the IPA asymptotic amplitudes for ionization of all the other occupied 
atomic states j: 

= ^° + E^%-., (1) 

with Aj^i the matrix element for a transition from the state i to the state j, caused by the 
outgoing electron, following photoionization of the state j. If j is a bound n's state (this 
will the most important case, since correlations with higher i states contribute beyond the 
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asymptotics [4]), the asymptotics of the IPA amphtude can be written as [10, 13] 



with the plane wave approximation for the wave function of the outgoing ejected electron, 
N^,g the normalization factor of the radial function of the n's electron ( ■0n's(^) — '4'n'si''") I \/47r; 
^n's — V'n's(O) , T = maZ , and 7 = — i(e ■ V)/m is the interaction operator between the photon 
with polarization vector e and the electron. 

The plane wave approximation 

V^f (r) = exp[i(P -r)], (3) 

with momentum of the outgoing electron P ^ r, is appropriate for Eq. (2) in velocity form. 
This corresponds to normalization of the continuum wave functions by the condition [11] 



J d^rjp*p{r)jPp>{r) = {2nfS{P - P') 



The amphtude is evaluated in Appendix A. 

If the photon energy well exceeds the binding energies of the bound states i and j, the 
matrix elements Aj^i can be represented as 

= ^^S,, , (4) 

with ^ = ma/ P {P the momentum of the outgoing electron) the Sommerfeld parameter of the 
final state interaction of the outgoing electron with the residual ion. The matrix elements 

Sj, = {j\\n{r - z)\i) , (5) 

(with z the projection of the coordinate vector r on the direction of the momentum of the 
outgoing electron), obtained in [1], describe the transfer of an electron from the states i to fill 
the hole in the state j of the positive ion with the hole in i state. One can write ln(r — z) — 
Inr + ln(l — t), with t = {P ■ r)/Pr and P the momentum of the outgoing electron. For 
states i and j with different angular momenta the matrix element (5) with Inr vanishes due to 
orthogonality of the angular parts of the wave functions of the states i and j. Thus when the 
states i and j have different angular momenta we can write 

% = (j|ln(l-i)|i). (6) 

Only such states contribute to the asymptotic amplitude, when i is not an s state (which we 
will assume), since the dominant asymptotic amplitudes Fj require j be an s state. Thus, we 
shall consider only the states j with quantum numbers £ — £z — 0, i. e. j = n\ 0, 0. Taking the 
direction of the outgoing electron momentum as the axis of quantization of angular momentum, 
we find that correlations occur only for i states with £^ = 0, since only these states are coupled 
by Eq. (6). 
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We may therefore write 

F, = F^ + iiT.,A,, , (7) 

where 

= FjSji, (8) 

with Sj^i given by Eq. (6) for photon energies, well exceeding the binding energies of i and j 
states (being suppressed otherwise [1, 6]). For example, as shown in [1], in ionization of the 2p 
electrons of neon by photons with energies of about 1 keV correlations with 2s electrons are 
important, while those with Is electrons are not; but by 10 keV correlations with Is electrons 
become important. We will write Aji as Aj omitting the index i. 

Now we restrict ourselves to the case of ionization of p states only, i. e. we consider the case 
i = n, 1, 0. We use the standard spectroscopic notation, e.g. the state with quantum numbers 
n, 0, is an ns state, etc. Then Eq. (6) can be evaluated as 

Sii = -^(C'.lV'^p), (9) 

with the first factor coming from the angular integration, ^0^^ and 0^/^ are respectively the 
radial wave functions in the field of the atom and of the ion with the hole in np state, and 

Note that since all the other electrons in initial and final states belong to different Hamil- 
tonians, there are nonzero overlap integrals between orbitals of different subshells. This makes 
the whole picture more complicated-see [12]. For example, in our case there are other channels 
for ionization of 2p state. In one of them the photon interacts directly with 2s electron, while 
the Is electron suffers shakeup into the hole in the 2s state of the ion. The photoelectron pushes 
the 2p electron into Is hole of the final state ion at the end of the story. The contribution of 
this channel to the total amplitude of ionization from the 2p state is thus -^2s^is,2p(02sl'^is)- 
|(02slV^L)l <^ 1 we can neglect this contribution with respect to the other terms on the right 
hand side of Eq.(l). 

We shall consider this very case, thus assuming that K^J^'sl'V^^^)] <^ 1 for all n ^ n' . Hence 
we find 1 — \{Vns\'^ns)\'^ ^ inclusive cross section with the sum over all 

possible states of the final ion coincides with the exclusive cross section in which the spectator 
electrons do not suffer transitions. In this approximation we must replace Eq.(9) by 

S,r = -^{Vn's\Vn,)- (10) 

Hence, 

An's = -^F^'Ai^l'sWnp) , (11) 

with the asymptotic IPA amplitude for ionization of an n's bound state. We shall omit 
the upper index IPA further. For correlations inside the same shell 



5 



For ionization of 2p states, i. e. n = 2, 

(V'LIV'U «^ -1, (12) 

for all atoms. This matrix element, calculated with Hartree-Fock wave functions, is —0.91 for 
Z — 5, with the value becoming closer to —1 for larger Z. The Coulomb value is — -\/3/2.Thus 
we can estimate that 

Ms = ^^2° ■ (13) 

For photon energies well exceeding the bounding energy of L shell, but not of the K shell, 
the correlation with the 2s electron dominates. Correlations with Is electrons are small at these 
energies [1, 4]. Correlations with other s electrons, if there are any, are small This happens 
for two reasons. The overlap matrix elements |('0n'slV'2p)l ^ l(V'2slV'2p)L 2, due to 

Eq.(12) and to the closure relation J2xs \ {i-^ls\'^2p) P = 1 (with summation over the states of 
both discrete and continuum spectra) Also the asymptotic IPA amplitudes F^,^ drop with n'. 
Hence at such energies the value (13) determines the scale of IPA breaking effects. 

However, at larger energies, greatly exceeding the binding energy of the Is electrons, Ais 
becomes comparable to A2s, and there is a large cancellation between contributions of the 
K shell and the other shells. Such cancellations were first found in the angular distributions 
of photoionization [3] and then observed at the amplitude level [5]. The calculations require 
knowledge of rather precise wave functions for the description of the bound electrons. 

For photon energies much greater than all the binding energies the asymptotic contribution 
of correlations to the amplitude is 

fa = ^-A. = Efi,A^, , (14) 

with the sum over all occupied hs states, where A^g is given by Eq. (11), for all occupied states 
j, where Aj is given by Eq. (8), namely 

A~. = (V'?|7|V'j)(0il lnr(l - tM) ^ (V'J|7|V'i)(V'j| lnr(l - tM). (15) 

There are only two active electrons in our analysis, while the others are just ignored. This 
corresponds to calculation of the amplitude for an inclusive process, since the sum over all 
possible final states of the spectator electrons in the final state ion, described in the initial 
state atom by the function ^ 

E(*/l*>' = 1 



3 Closure condition and sum rules 

Using closure we can write simple relations for amplitudes like in Eq. (14), but now summed 
over various states j. We will do it in a 3d formalism. We will sketch these relations here and 
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complete their proof subsequently. The closure condition for the wave functions of the initial 
state nonrelativistic Hamiltonian can be written as 

E = 1 = E + E IMiM + E IV'.c)(V'.el, (16) 

J j J 3c 

with j and j* labelling the occupied and unoccupied states of discrete spectrum correspondingly, 
while jc are suitably normalized continuum states. The closure condition can be represented as 

/7|^^QW^g(0+ E <wWn'Am'(r') = ^(r-r'), 

V n',£',m' 

for \1'q which are asymptotically plane waves as in Eq.(3). For any state j, occupied or not, we 
can write, generalizing Eq. (15), 

lit fit 



Wc may evaluate Eqs. (14), (15) using Eqs. (16) for the sum J2j \'4'j){'4^j\ over occupied bound 
states. 

fd = T-T*-n, (18) 



with 

T= EjAj = (^°|7lnr(l-t)|V'.) = lnr(l - t)|^,) (19) 

the sum over the complete set of states j, while Td and are the sums over occupied and un- 
occupied states j of the discrete spectrum correspondingly, with Tc the sum over the continuum 

states. 

We can write Eq. (18) in the form 

T = Td + n, (20) 

with 

T, = f, + T*, (21) 

the sum over occupied and unoccupied states of the discrete spectrum. 

We may now separate the continuum amphtude Tc defined below Eq.(19), for which the 
continuum states jc may be labeled by asymptotic momentum Q, in two parts 

Tc = Tci + Tc2] T = Td + Tci+Tc2. (22) 

Tci will sum states for which q <^ P, Tc2 states with Q P. More precisely, we pick a Qq for 
which p S> Qo ^ Tc, with Tc the characteristic momentum of the bound state (one can assume 
Tc ~ maZ), and define as the sum over states Q < Qo, Tc2 as the sum over Q > Qq. 

We may show (see Appendix B) that for states in Tc2 the wave function ipq may be replaced 
by a plane wave: ipq = exp [i(Q ■ r)]. We may also show (see Sec. 4) that the sum over plane 
wave states Y,q>Qo \'4^q){'4'q\ iii 2^c2 niay be extended to a sum over all plane wave states J2q, 



that is, the sum over plane wave states, Q < Qq makes no contribution in Tc2. But J2q IV'q) H^qI 
over plane wave states is a sum over a complete set of states (=1), and therefore for Tc2 one 
again obtains Eq. (19), i.e. in the asymptotics 

T ^ T,2. (23) 

In Appendix C we show that the amplitude T'c2 has the same asymptotics as T^.. Hence, 
while from Eqs. (19), (20) and (21) 

T = fa + T* + T,^+T,2, (24) 

Eq. (23) implies 

o = f, + T; + Tei; frf = -T;-r,i, (25) 

giving an alternative way to calculate which, as we shall see, has some advantages. (In fact 
all these are asymptotic amplitudes, and so for Eq. (25) to follow from (23) and (24) we must 
show that these are all amphtudes of the same order. We shall do this in Sec. 4.) In Sec. 4 we 
shall make exphcit calculations of the amphtudes of Eq. (25), determining from and Tci, 
as well as further discussing its direct calculations in Sec. 5. 



4 Formalism for particular amplitudes 



a In 3d formalism all of the particular amplitudes of the previous section are of the form of 
Eq. (14), except that the summation J2j in is replaced by summation in for unoccupied 
bound states, and by integration / in Tc for continuum states with Q < Qo and Q > Qo 
for Tci and Tc2 correspondingly. Also, T, corresponding to summing over the complete set of 
states, was given explicitly in Eq. (19). 

We can also represent the amphtudes in terms of radial functions. Using Eqs. (2) and (10) 
we can write 

Td = E^l^K'sWn'sl^np) (26) 

rn r- ^, 

and correspondingly for partitions and summing over occupied or unoccupied bound 
states n's. 

Continuum radial wave functions for s states (as well as those for a nonzero value of €) are 
normalized by condition 

with V''P,s(^) = ^■''p,s(^)/(2'?'"p)^^^, 'Vpiy) — V'p^s (?') /2p+terms with nonzero values of I. The 
closure relation is 

/ dei,Z[r)i^lXr') + E C*.(r)V'n',.(r') = b{r - r'). 

n',s 
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Thus the calculation of Td proceeds in the same way again (see Appendix A), and yields 
Eq. (26), except that for the radial functions we write de instead of Xln'; ^-^d replace N^,g 
by N^g. (The integration over energies has been extended to infinity, with the energies exceeding 
strongly the binding energy of the ionized state, including those for which e > Ql/2m, providing 
a negligible contribution.) 

For Tc2, where Q > Qo, we may argue (Appendix B) that in Fq, corresponding to Eq. (2), 

■0Q is to be replaced by a plane wave. Thus = ^^^S{P - Q). Hence, with plane wave for 
i/jQ we may extend the integration over Q to include Q < Qq <^ P ^ since there is no contribution 
from this region. On the other hand, we can evaluate explicitly (see Appendix C) 



" m P"^ ^ ' 

As we have already seen, T = Tc2- Note that we have also now shown that all these terms are 
of the same asymptotic order, which we had needed to prove Eq. (25). 

Note that we can write Eq. (26) in the form 

Zl-^i2'S'n's,np + j J^^^Q^Sq^np = 0, (28) 

with Sj^i defined by Eq. (5), F^?] arc the high energy IPA photoionization amplitudes [13] while 
Fq^ is the brcmsstrahhmg amplitude in the tip region [14]. Following the previous analysis, the 
integral in the second term involves all values of Q. However it is saturated by Q ~ r = maZ. 

We can now write Eq. (25), expanding the radial function V'^p(r) in terms of the functions 
Vxsir) {x = n',e): 

■oc 

^np(^) = Y.^n'sir)an's,np + J deiljlg{r)a^s,np : (29) 



with 

an's,np = ii^n'sl^np) ' (^es,np = (^n'J^np) ' (30) 

while closure can be written as 



oo 



n' 

(Note that axs,np = —'^/V^Sxs,np, with Sxs,np given by Eqs. (6) and (9)). Note that integrals 
on the RHS of Eqs. (29) and (31) are saturated by energies of the order of the np electron 
binding energy. 

The ratio of correlations of the np state, with n's and es, to correlation with the ns state 
may be described by the factors 

_ Nn's^n's,np _ -^gs^e.n ,ncy\ 

•i'n's,np ~ jvfr ^ ' •^es,np ~ Ajr ^ ' ^ ' 

^^ns^'ns,np ns^n.n 
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We can write Eq. (28) in the form 

oo 

^n's,np "I" J d£Xss,np — . (33) 
n' 

We must calculate the physical value of the total correlation, relative to the correlation 
between np and ns electrons, i.e. the sum of relative correlations over occupied states n, 

^ ! -^hSjUp 1 (34) 



X. 



ph 



n 



measuring the total amount of correlation relative to the intrashell ns, np correlation. 

5 Calculation of the physical value Xph 

Using Eq. (32) for Xn's,np one can see that there is a tendency of cancellation of correlation 
effects for the 2p electrons. Employing Eqs. (12) and using Eq. (31) we find |a„'s2p| ^ 1 for 
n' ^ 2. Since the normalization factors N^,^ drop when n' increases, the contribution of these 
states to the RHS of Eq. (34) can be neglected for n' > 2. The case of n' = 1 requires additional 
analysis since A^is > A^2s- One can see that the contributions of n' = 1 and n' = 2 to Xph have 
different signs; from Eq.(12) a2s,2p < 0, while ais,2p > 0, since the function has no nodes, 
while the function i/j2p is non-negative. Hence, there is a partial cancellation between these 
terms. We shall discuss the extent of this cancellation subsequently. 

Now we try to calculate Xp^ in another way, i.e. as the negative of the total contribution of 
low energy unoccupied states, bound (xd) and continuum {xc): 

•^ph — ^ y -^n^n — -^d -^o -^d — ^ ^ •^n*s,np i -^c — / d£Xgs,np (35) 

h n* 

with n* labelling the unoccupied states of the discrete spectrum. 
5.1 Discrete states 

It is known that for n' ^ n the dominant region of coordinate space in the integral {'i'n's\''Pnp) 
determined by the characteristic size of the np state. These values of r are much smaller than 
the characteristic size of the n's state. The dependence on the energy of the n'-ih state in the 
Schrodinger equation for ip^^,^ can be dropped [10], and the only n' dependence of this matrix 
element is contained in the normalization factors Nn's- Thus, from Eq. (32), the ratio 

Xn's,np — ' i^^) 



ns 



To estimate A^^*^ for the unoccupied states with principal quantum numbers n* we use the 
quantum defect approach, in which the binding energy E^' of the n's state is 

^ " 2(n' - A„0^ ' ^^^^ 
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with V = ma"^, and 

An' —>■ A as n' — > oo . 
Here the quantum defect A does not depend on n'. 

Combining Eq. (37) with the Fermi-Segre formula [15], 



(38) 



we obtain 



with 



n- An 
n' - An 

l-AL 



(39) 

(40) 
(41) 



We neglect the derivatives A^,, as justified below. Applying Eq. (39) to the highest occupied 
state fih of known binding energy, we find its quantum defect Ah. Alternatively, for the lowest 
unoccupied level of s electrons, n* = n^^ + l, we identify the quantum defect A„ with the limiting 
value A, defined by Eq. (38). The latter can be extracted from the results of [16], on the phase 
shifts with respect to Coulomb values S{E), since A = 5(0)/7r. The values of A/j, calculated 
by using the normalization factors N^/g given in [17], and A, are shown in Table 1. We will 
assume that A„ for all unoccupied states is independent of n, An — A. The comparison of Ah 
and A indicates the type of error that is being made. Relatively small value of the difference 
A/j — A justifies neglect of the derivatives A^,. 



z 


A, 


A 


5 


0.96 


0.76 


7 


1.23 


0.95 


10 


1.44 


1.27 


14 


2.00 


1.69 


18 


2.31 


2.04 


32 


3.04 


2.74 


36 


3.28 


3.06 


50 


3.96 


3.37 



Table 1: The quantum defects of the highest occupied bound states Ah, and the asymptotic 
values A as defined by Eq. (38), obtained from [16]; Z is the nuclear charge. 



The total contributions of unoccupied discrete levels can be written as as 



(42) 
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with Hh the principal quantum number of the highest occupied state. The summation over 
unoccupied states in Eq. (42) can be carried out by using the formula [19], 



oo 1 -1 rifi 

with ip{a) = r'(a)/r(a), where r(a) is the Euler gamma function, fih + 1 is the principle 
quantum number of the lowest unoccupied state. 



5.2 Continuum states 

Equations (37) and (39) reflect the Coulomb-like behavior of the excited states of the discrete 
spectrum at n' — > oo. Thus for the continuum states with £ = we can write 

XOs,np = 7^ • liin in'^Xn's,np), (44) 
ZIq n'— >oo 

with lo = ma'^/2 = 13.6 eV, or 

AZIo 

XOs,np - ■ (45) 



ns 



Since the only characteristic energy is the binding energy Enp < of the ionized state, we 
may suppose that the integral in the last equality of Eq. (35) is saturated by £ ~ \Enp\- If we 
suppose that Xss,np = Xos,np for s < \Enp\ Xss,np = for £ > \Enp\ we find 



AZIqIE, 



np\ 



(46) 



This is clearly a fairly crude estimate, but we will see that it is consistent with results from 
direct calculations of Xph. 



5.3 Results for Xph 

Note that in the quantum defect approach > 0. The assumption (46) provides also x^ > 0. 
Hence the values Xph arc negative. This means that the total correlation effect in the amplitude 
has a sign, which is opposite to that of correlation inside the same shell. Employing Eqs. (42), 
(43) and (46) we find the values of Xph shown in Table 2. 

Our results are in good agreement with the results x"^^^ obtained by direct summation of the 
correlations with occupied shells in the photoionization amplitude [5]. Results of [5] obtained 
by inclusion of correlations with 2s and Is shells in nitrogen and neon arc Xph = —0.18 and 
Xph = —0.11 respectively. Since the results of [5] are in good agreement with those of [4] for 
angular distributions, our results agree with those of [4] as well. For ionization of the 3p state 
in argon, "shell by shell" calculation [5] gives Xph — —0.14, also in agreement with the result of 
the present work. 
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Zj 


n 


Xd 






a; 


r 



2 


(j.(j8 


A 1 O 


A OA 




( 


o 
/ 


n nix 


u.uy 


-U.i4 


-U.18 


iU 


2 


U.U4 


A A'7' 

U.U7 


All 

-U.ll 


nil 

-U.ll 


14 


2 


U.Uo 


A 1 

U.lo 


A O O 

-U.22 




1 A 

14 


3 


A AO 

U.U8 


A 1 O 

U.12 


A OA 

-U.2U 




1 O 

18 


o 


A C\P 

U.Uo 


A OA 

U.2U 


A O/^ 

-U.2o 




1 8 

lO 


Q 

o 


U.UO 


U. lU 


-U. lO 


n izL 

-U . 1^ 


32 


3 


U.Ul 


U.18 


-u.19 




32 


4 


U.13 


U.ll 


-U.24 




36 


3 


U.Ul 


U.19 


-U.2U 




36 


4 


U.U8 


U.ll 


-U.19 




50 


5 


0.06 


0.13 


-0.19 





Table 2: The values of Xc and Xp/j as defined by Eq. (35), where n is the principal quantum 
number of the ionized np state, Z stands for nuclear charge. The values x**" presented in the 
last column are the results of direct summation over occupied states - Eq. (34). 



Now we estimate the total contribution of correlations to the amplitudes and cross sections 
of photoionization of p states. We use the estimate {ipnp\ipns) ~ —1 for all n (see Eq. (12)). Pre- 
senting the ratio of IPA amplitudes in terms of normahzation factors -Fn,o,o/-^n,i,o = ^l/V^^p: 
we find 

Fn,i,o = <i,o(^l + ^^p.J, (47) 
and thus the cross section for ionization of np state beyond IPA is 

/ iV \ 

^np = (^np[^ + ^^phj, (48) 

with (T°p standing for the IPA values. 

Using the numerical values of the normalization factors [17] we find that the total correla- 
tions diminish the values of cross sections of photoionization of 2p states in nitrogen and neon 
only by about 2.5%. In contrast, inclusion of correlations only with the 2s shell would increase 
the cross sections by 18% and 22% correspondingly. The full cross section for ionization of 3p 
states in argon becomes smaller by 1.8%, while it becomes larger by 12% if only correlation 
with the 3s shell is considered. 
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6 Calculation for hydrogenlike functions 

and limitations for the many-electron atoms 



We can also make explicit calculations of correlations in the case of Coulomb wave functions, 
and also in using an effective charge approximation for screening. All results for Coulomb 
functions can be obtained analytically. The results for Xn's,np do not depend on the values of 
nuclear charge Z. 

Starting with ionization of 2p electrons we obtain the values of the parameters that are 
presented in Table 3 (upper indices C indicate that the quantities are calculated in the Coulomb 
field of the nucleus) . 



n' 


c 


6 


1 


-1.58 


-1.26 


2 


1.00 


-0.02 


3 


0.041 


1.00 


4 


0.015 


0.04 



Table 3: Parameters for ionization of 2p and 3p states, obtained by using Coulomb functions. 



One can see that in this case the correlation with the Is state is about 50% larger than that 
with the 2s state. For the contribution of the continuum with s <^ E one can obtain by direct 
calculation 

^-2P = 77T7^; $(£)=exp(-2ei(arctan(2/ei))-2/ei)), ^ r/e ^ ^fl^e; $(0) = 1, 



where 



(49) 

C = 2 hm n'^x^,^^,^ - 0.78 (50) 

n'— »oo ' 

is obtained by using the well known Coulomb wave function for the bound n's state for n' ^ 1 
[10]. This gives Xc — 0.52, in agreement with Eq. (33). 

Consider now ionization of 3p electrons - see again Table 3. The correlation with Is and 
3s states cancel to a larger extent than in the case of the 2s state. The contribution of the 
continuum is now Xc = 0.22. 

Now wc analyze the situation for more realistic atomic models. Wc shall compare the values 
Xph calculated is the unscreened and screened Coulomb fields. Here we obtain Xph as the sum 
over occupied states, taking into account that, as we showed above, for ionization of a np state 
only correlations with ns and Is states are important. Using Eq.(32) we can write in this 
approximation 

with Xis^np < 0, as shown above. As one can see from Table 3, the Coulomb values > 1. 

Now we show that for the screened Coulomb values < and thus the Coulomb 

values x'-^ can be used as the lower limits for the physical values Xph. 
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We can calculate the screening effects, assuming that the initial electrons are described by 
the Coulomb functions with effective values of the nuclear charge Zni = Z — [10, 18]. In 
this approach we find for ionization of 2p electrons 

^i.,np - v-[z,J \2Zu + zJ V 2Z ) 3^2. ^^^^ 

with 77 = 1 if screening is neglected and thus Zis — Z2S — Z^^ — Z. In the lowest order of 
expansion in powers of bni Eq. (52) provides = 1 + with 5 = — 5is/3 — 552^/3 + 2^23- 
If a small influence of the electrons in the higher states on the values of 8ni is neglected, 8 is 
the same for all atoms with the totally occupied K and L shells. Using the values 8\s — 0.35, 
= 3.25 and = 4.75, [18], we find 5 < 0. Thus r] < 1, and indeed 

Xphys > X -I (^"^) 

while x*" < 0. Using Eq.(52) for neon with these values of we find 77 = 0.676, providing 
Xis,2p = —1.068 and thus Xphys = —0.068, with |a;p/jj,s| smaller then that shown in Table 2. Note 
however that this value is a result of subtraction of two much larger values. Putting tj = 0.702, 
i.e. increasing it by 4% we would find Xp/jy^ = —0.11, in agreement with the data in Table 2. 



7 Exclusive and inclusive processes in neutral atoms and 
ions 

Instead of photoionization of neutral atoms we can consider photoionization of ions. Since the 
cancellation of correlations is due mainly to cancellation between correlations with Is and with 
ns electrons, a hole in either the Is or ns shell breaks this balance, and the net correlations 
will greatly increase. This was observed earlier [3] for the cases of nitrogen and neon. A hole 
in other shells will not influence strongly the total correlation. 

One is often interested in inclusive and exclusive photoionization, as in ionization of a 2p 
state, but perhaps also exciting other electrons. The theory of such processes was much studied 
[20]-[22] in the case of shake off and shake up. 

As already noted, our discussion considering only two active electrons, was inclusive in its 
treatment of the spectator electrons, it could be exclusive if overlap integrals between initial 
and final spectator states were considered. Note that in principle one should have also, in the 
presence of correlations, include other electrons as active, capable to undergo further excitation 
or ionization beyond the shakeoff/shakeup mechanism. We can try to estimate the magnitude 
of these various mechanisms. 

Assume that the bound electrons, moving in a certain self-consistent field , find themselves 
in another field after the electron is ejected. Photoionization of an n, i state can be followed by a 
transition of an electron n', £' to the state n*,£*. One should consider the process simultaneously 
with ionization of n', £' state followed by a transition of an electron n, i to the state n*,£*. The 
asymptotic IPA amplitude of the process, without correlations, is 
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with ip and the functions in the fields of the atom and of the ion with a hole in n, ^ and n', 
states for the two terms in the right hand sides (RHS) of Eq.(54) correspondingly, is the 
asymptotic IPA amplitude for ionization of nl state. One needs t — I' , or i* — i otherwise the 
matrix element vanishes due to orthogonality of the angular parts of the wave functions. 

Correlations provide another mechanism of the process in which n,£ electron is ionized 

by direct interaction with the photon, and in the next step the photoelectron excites the n'i' 
electron to n*i* state. One should include possible permutation of the n£ and n'£' states. The 
amplitude, which includes the correlations can be written as 

= + F^e^n*e*,n'e' — Fn'£'-^n*e*,ne, (55) 

with J-'^ given by Eq.(54). The correlations can cause these transitions even if £* coincides 
neither with £ nor with £', and the shakeup mechanism can not contribute. 

Note that the contribution of correlations on the RHS of Eq.(55) is written omitting the 
terms containing as additional factors the overlap matrix elements of the type {(t>n'e\'4^n'e) ■ We 
neglected such terms in particular calculations through the paper-see Sec. 3. Inclusion of such 
terms would not alter the asymptotic energy dependence of the amplitude. Hence we shall use 
Eq.(55). 

Consider, for example, photoionization of Be with the final state ion containing electron 
excited into 2p state. In this case n = l,n' = n* = 2,£ = £' = 0,£* = 1. Both terms on the 
RHS of Eq.(55) turn to zero, providing = 0, and thus 

^=F°A2p,2.-F°A2p,i,. (56) 

Now we study the the relative role of the shakeup and correlation mechanisms of the process, 
for various relations between £, £' and £*. To obtain the energy dependence of the contributions 
on the RHS of Eq.(55) one can employ that in the asymptotics 

F^, ^ cv-^'+^y' An*.*,n. ~ (57) 

The estimation for is well known [10]. The estimation for An*e*jd is the consequence of 
Eq.(4), in which the matrix element of A between the bound states is proportional to the factor 

As we have seen for £ ^ £' ^ £* only correlations contribute. Turn now to other cases. If 
£^ £' — £* the second term on the RHS of Eq. (54) vanishes and the IPA amplitude is 

Using Eq.(57) we find that the second term on the RHS of Eq.(55) drops faster than and 
hence the asymptotics is 

^ — F^e{<Pn*e'\i^n'e') — F^'i'An*i',ni- (58) 

Further analysis depends on relation between £ and £'. 

For £' = £*<£— l(for example, n's n*s and n'p n*p transitions in ionization of d 
states), the correlations determine the asymptotics of the process since the first term on the 
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RHS of Eq. (58) drops with energy faster then the second term. For = = (for 
example, n's — > n*s transitions in ionization of p states) states the two terms terms on the 
RHS of Eq.(58) behave with energy in the same way. However the first term is proportional to 
the overlap matrix element, which is usually small. If it is the case, the correlations dominate 
the process. 

Similar analysis shows that for C' = C* > i — 1, including the case i = i' = i* (for example, 
n's — > n*s transitions in ionization of s states) asymptotics is determined by the shakeup 
mechanism, described by the first term on the RHS of Eq.(58). However, at finite energies, 
where the experimental data is available, interplay of the shakeup and correlation terms appears 
to be important [9], [23]. 

To obtain the cross section for the inclusive process one should sum the squared amplitude 
given by Eq.(55) over n' and n* , depending on the conditions of experiment. 

8 Summary 

We have calculated the IPA breaking correlation corrections to the high energy photoionization 
amplitude, focusing on ionization of p states. Instead of carrying out summation over occupied 
states, we employed the closure results for summation over all states of the spectrum. We 
showed that the sum over all states of the spectrum is equal to the contribution of its high 
energy part. Therefore there is total cancellation between contributions of discrete and low 
energy part of continuum spectra. This provided identities, involving the asymptotics of the 
amplitudes of photoionization and of bremsstrahlung amplitudes at the tip region - Eqs. (28), 
(33). 

We calculated the sum of correlations with the occupied bound states as the negative of 
the sum over the unoccupied bound states and the low energy continuum states. We made 
conclusions in a simple model, based on the general features of the bound states with large 
principal quantum numbers n. In this approach the sum of the occupied state correlations has 
a sign, which is opposite to that of the correlation inside the same shell. In spite of a crude 
model for the continuum, the results are in good agreement with those obtained earlier in direct 
calculation (Table 2). 

We have shown that there is a general tendency of cancellation for the correlation effects. 
We demonstrated also that calculations with Coulomb functions give limits for the correlation 
effects in screened atoms, Eq.(53). 

We showed also that the correlations beyond shakeoff and shakeup effects are important 
in inclusive processes, where photoionization is accompanied by excitation of other electrons. 
The relative role of the shakeup and correlation mechanisms was found to depend on relations 
between orbital momenta i and £' of the removed electrons and orbital momentum i* of the 
excited electron. In some of the cases the correlations dominate in the process. In the particular 
cases, for which experimental data are available , interplay of the two mechanisms is important. 
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Appendix A 

We calculate the asymptotic amplitude (2) for the bound n's states, following the approach of 
[13]. Note that it is determined at small r ~ 1/P. Thus we can use expansion of the function 

ilJnsir) at r ^ 0, 

V'n.(r) = ^(l + a,r)exp(-Ar); A > 0, A ^ 0. (A.l) 

v47r 



Here the last factor has been introduced to insure the convergence of the integral in the in- 
termediate steps. We have kept two terms of expansion in powers of r in brackets, since the 
lowest one, as we shall see below, vanishes at A = 0. (Higher terms in r would contribute in 
higher terms in 1/P). The parameter on the RHS of Eq. (A.l) should is equal to the first 
derivative of the function '0ns (r), as determined by the first Kato cusp condition [24], being 
as — —maZ = — r. Since 



/ 



and re ^'^ — (— ^) e we obtain Eq. (2). 



Evaluation of the matrix element in Eq.(2) corresponding to a continuum state j with 
asymptotic momentum Q <^ P can be done in the same way, with the same form of expansion 
of %l)Qs. Since only s waves contribute, we can write 

where e is the energy of the continuum electron, 

47rV2 ]\jr 

Xes = p, " , (A.4) 

with = ipl^^O), where the upper index r again denotes the radial part of the function, 
yielding Eq.(2). 



Appendix B 

Now we evaluate the amplitude Tc2 defined, following Eq. (17), as 

(eP) r d^Q 



Tc2 = 



m 
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with integration over Q ^ t, and with 



X{Q) = {^Q\lnr{l-tM). (B.2) 
Using Eq. (9.6) of [10], and taking the first iteration, wc obtain 

irfli^Q) = {27rf6{P-Q)+h{Q) (B.3) 

with 

The first term on the right hand side of Eq. (B.3) immediately gives 

(eP) 

Tc2 = ^X(P) , (B.5) 
m 

leading to Eq. (23) due to Eq. (19). 

Now we show that the second term on the right hand side of Eq. (B.3) provides higher order 
terms of expansion in powers of P~^. 

We put 2mV{q) = Tcv{q^), with Tc being of the order r <^ P. One can see immediately that 
the regions \Q — P\ ~ P lead to corrections of the order 1/P. The vicinity of the point Q — P 
requires special analysis. Near this point we use the well known relation 

^ = P.V.--m5{0). (B.6) 



X + IU X 

For the first term on the right hand side of Eq. (B6) the result of integration over the angles leads 
to a function of (P — QY, i.e. to an even function ol P — Q. Together with the denominator 
P — Q this leads to an odd function of P — Q in the integrand of the integral over Q, providing 
contribution of the order ~ TcX{p)/P. Contribution of the same order comes from the whole 
interval Q ^ P. In a similar way one can see that the second term on the right hand side of 
Eq.(B6) also contributes only beyond the asymptotics. Thus, indeed we can neglect the second 
term on the right hand side of (B.3). 



Appendix C 

In order to calculate we must evaluate the matrix element defined by Eq.(35). It is expressed 
also by Eq.(30 of appendix B. Since Q » r we first describe '4>q by a plane wave (3). Then 

X{Q) ^ J dVe-'^nn{r-z)iJi{r). (C.l) 

For p states we can write 

i^iir) = /J^V'I(r) (C.2) 
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(recall that we need only the states with 1^ = 0). Since we shall need ipi{r) at r ~ l/Q, we 
can put iplir) = Nlre'^'' (A ^ 0) in Eq. (A.4), and thus ^i{r) = ^3/47r iVfre-^^ We use 
^g-iQz _ iJjge-'i'Q^ (with Dp = d/dp) and use the parabohc coordinates 

^^r + z, T) — r — z, (f) — axctsin{x / y) , 

so that 

2 ' ^~ 2 ' 

and 

dV — — d^dr]d(j) ■ 

Thus we find 

X{Q) = ^^^^^^ i^cji?. Jg (C.3) 

with 

with 7s f» 0.578 the Euler constant, leading to asymptotics 

Note that the asymptotics Q"*^ of the function X{Q) is due to the logarithmic term in the 
integrand. Replacing it by a constant, we would obtain the asymptotics Q^^. Indeed, replacing 
ln?7 by a constant, we would immediately find X{Q) = due to the operator Dx. Thus we 
must include the next term of expansion of the function ipi{r) in powers of r in a way similar to 
Appendix A. This gives X(Q) ~ DqDI ■ 1/(X'^ + Q^) = It can be shown that the higher 

order corrections to the function x/jq contribute only to higher orders of expansion in powers of 

1/g. 
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